syms x y a b c d;
[n1,d1] = numden(sym(tan(3/5))) 
[n2,d2] = numden(x/y + y/x)
A = [a, 1/b;1/c d];
[n3,d3] = numden(A)
